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Abstract 

We review both the counting rule and the influence of the evolution in Q 2 for the 
large XBj behaviour of the valance quark distribution functions. Based on a factor- 
ization procedure we present a more general perturbative treatment to compute this 
behaviour. A complete analysis is performed in the scalar (^-theory for the parton 
distribution function of the "meson", which shows that logarithmical corrections arise 
from the distribution amplitude and that the reference momentum square Qq is fixed 
by x B j- 



1 Supported in part by Deutschen Akademischen Austauschdienst and by Department of Energy contract 
DE-AC03-76SF00515. 



1 Introduction 



Deep inelastic lepton-hadron scattering (DIS) experiments have shown that hadrons consist 
of point-like particles called partons. These partons were identified as quarks and gluons, 
which are the fundamental degrees of freedom of a relativistic SU C (3) gauge field theory, 
quantum chromodynamics (QCD). For instance, corresponding to this quark-parton model 
a nucleon is formed by a colour singlet state that consists of three valence quarks as well as 
sea quarks and gluons. 

Since the bound state problem remains unsolved at present it is not possible to compute 
the DIS cross section completely from the first principle of QCD. However, with the help of 
the operator product expansion |]J it could be shown that the hadronic part factorises as a 
convolution of perturbative computable Wilson coefficients and nonperturbative parton dis- 
tribution functions 0].R These distribution functions are defined as forward matrix elements 
of string light-cone (LC) operators ||, J|]. Moreover, from the renormalization group (RG) 
equation of these operators one obtains the Gribov-Lipatov-Altarelli-Parisi (GLAP) equa- 
tion ||, which controls the evolution of the distribution functions in Q. However, it remains 
a theoretical non-predictable input, the distribution functions at a reference momentum 
Qo ~ 1 GeV. 

At large x B j the theoretical situation is much better. Therefore, it is possible, with the 
help of the naive quark-parton model, to derive a counting rule || that predicts the power 
behaviour of the distribution functions in (1 — x B j) for both mesons and baryons |7j]. These 
predictions will be reviewed in Section 2. 

In the Section 3 we define the distribution functions as expectation values of LC operators 
and discuss, with the help of their RG equation (or GLAP equation), the change of the large 
xbj behaviour by evolution in Q. 

In Section 4 the distribution functions will be factorized as a convolution of distribution 
amplitudes and a perturbative computable hard amplitude, which is computed for the "me- 
son" distribution function of the scalar 0^ -theory in six dimensions. Apart from the leading 
power behaviour the convolution also provides logarithmical correction in (1 — x B j) that 
arises from the evolution of the distribution amplitude, and also the value of the reference 
momentum square Qo as a function of x B j. 



2 Counting rule 

It is well known that the hadronic part of the unpolarized charged lepton DIS is described 
by the hadronic tensor 



W^(x Bj ,Q 2 ) = (-g^ + ^ Fl (x Bj ,Q 2 ) 



+ {p, - (p* - ^r) f*(*bi, Q 2 )/P q , (i; 



2 At first this factorization was performed in terms of the moments (expectation values of local operators) 
from which the mentioned convolution can be obtained with the help of the Mellin transformation. 
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where P is the four-momentum of the hadron, and q denotes the four-momentum transfer. 
Here Fi(xBj,Q 2 ) and F 2 (xBj,Q 2 ) are the structure functions which depend on the Bjorken 
variable < x B j = —q 2 /2qP < 1 and the momentum transfer Q 2 = —q 2 . 

In the Bjorken region, i.e., large momentum transfer Q and fixed XBj, the two structure 
functions can be expressed according to the quark-parton model by the non-perturbative 
quark distribution functions qi(xBj,Q 2 ) and anti-quark distribution functions: (ji(xBj,Q 2 ), 

Fx(x Bj ,Q 2 ) = F 2 (x Bj ,Q 2 )/(2x B] ) = \Y.^ (k(xb 3 ,Q 2 ) + q l (x B] ,Q 2 )) , (2) 

i 

where i — u,d, ... is the flavor index, and e% is the corresponding electrical charge. 

As mentioned, at large xsj it is possible to compute the parton distribution functions. 
But one deals with two large scales, Q — > oo and 1/(1 — x B j) — * oo, so one has to clarify in 
which way the two limits will be performed. With the help of the hadronic final state mass 
Ml = (P + q) 2 = Mjj + (1 — x B j)Q 2 /x B j > M%, where M H is the mass of the hadron, one 
finds the following classification: 

1. (1 - x Bj )Q 2 -> 0, M 2 X = M 2 H , 

2. (1 - x Bj )Q 2 -> finite, Ml = M 2 H + finite, 

3. (1 - x B j)Q 2 -> oo, M 2 X -»• oo. 

In this paper we consider only the third case, for which the factorization in Q via the operator 
product expansion remains true so that, as usual, higher order twist terms can be neglected. 

Now let us review the counting rule for the large xsj behaviour of the structure function. 
Corresponding to the naive quark-parton model we have the struck quark with momentum 
p, and the spectator quarks carry the momentum P — p: 




M 2 + k 2 l r \ „ /,„ M 2 + k\ 

1 - xb,-)P, 



V = ( x Bj P+, „ x , k± ] , P- P = |(1- x Bj )P+, —^r, -k ± ] . (3) 



Here we have introduced light-cone coordinates, i.e., p = (p + = p +p 3 ,p_ = p ~P3,p±), so 
that XBj = P+/P+. The essential assumption for deriving the counting rule is the requirement 
that the spectator quarks have a non- vanished finite mass M.. Then from the kinematical 
identity for the hadron mass Mjj = ((P — p) + p) 2 it follows that for xbj — ► 1 the mass of 
the struck quark 

2 2 x B] (l-x Bj )M 2 -x B] M 2 -k\ M 2 + kj 
M =p = ~ — (4) 

1 - X B j 1 - X B j 

is far off shell. Note that this behaviour also remains in the infinite momentum frame. 

The power behaviour of the structure functions at xsj — » 1 follows directly from the far off 
shell behaviour of the wave function. The iteration of the Bethe-Salpeter equation provides 
that this behaviour is determined by the minimal number of gluon exchanges required to 
stop the spectator (see Fig. 1 for a nucleon). Thus, without taking into account the spin of 
the partons one gets the following counting rule J7[: 

F 2 (x Bj ) ~ q{x B j) ~ (1 - XBj) 2113 ' 1 for x Bj -> 1, (5) 
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where n s is the number of spectators. 

Corresponding to the counting rule (J|) a nucleon has a power behaviour of 3, which is 
consistent with the experimental results. A meson should have a power behaviour of 1. But 
it was shown in Ref. || that spin effects provide an additional suppression so that the power 
is 2. It is interesting to note that experimentally one finds for the pion the power behaviour 
is closer to 1 than to 2 [jj; so that the situation about this prediction is not quite clear (see 
also Ref. § ). 



3 Improving the counting rule result by evolution 

The predicted behaviour of the structure function at large xsj will be modified by the evolu- 
tion of the distribution functions and by perturbative corrections of the Wilson coefficients. 
Although at Xsj — » 1 higher order corrections are important ||10|| , in this paper we discuss 
only the modification by the evolution at leading order. Interesting results from pertur- 
bative all-order investigations can be found in the literature for both the evolution of the 
distribution functions [11] and corrections of the Wilson coefficients [|12 . 



First let us sketch why the distribution functions are defined as expectation values of LC 
operators and why the GLAP equation arises from the RG equation of these operators. The 
dispersion relation 

W^{x Bj ,Q 2 ) = ^-lmT^(x Bj ,Q 2 ), 

47T 

T^xb^Q 2 ) = tjd 4 x< P\T{j„(x)j„(0)}\P > e***, (6) 

provides that the hadronic tensor can be expressed by the imaginary part of the forward 
compton scattering amplitude T^ u (xBj,Q 2 )- Furthermore, the leading Q 2 behaviour of 
W^ u (xBj,Q 2 ) is determined by the light-cone singularities of the time ordered product 
T{j fi (x)j u (0)}. Thus, the non-local light-cone expansion of this operator product |13| pro- 



vides a definition of the quark distribution functions in terms of gauge invariant string 
operators that are a resummation of local leading twist operators (see, e.g., flip. In the 
axial gauge xA = these operators are simplified to be bi-local: 



g,(£,Q 2 ) = e 2 ^ x <P\U-K~x)xiA{™)\P>\^ =Q 2, 0<z<l, (7) 

g,(x,Q 2 ) = -J |V 2 «(^ < P|^(-«z)57^ f («z)|P >^ =Q 2, < x < 1. (8) 



Here the more arbitrary light-like vector x arises from the projection of x onto the light-cone 
and can be chosen as x = (x + = 0, x_ = 2, 0j_) so that xP = P + . 

In the following we consider only the evolution in the flavor non-singlet sector. For 
technical simplification we define formally the flavor non-singlet operators 

O a (x,fi 2 ) = J ^e 2 ^ x U-^lKM^)^ ( 9 ) 
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where A"- is a generator of the flavor group. If we keep in mind the forward case then the 
renormalization group equation can be formally written as 

^O a (x, /i 2 ) = f l ^P(x/y; a s ^ 2 ))O a (y, /i 2 ), (10) 
a/r Jx y 

where the integral kernel is known as perturbative expansion in a s , 

P(x; a.) = C F ^- { ,] +X l + h(l - x)) + 0(a 2 ), where C F = 4/3. (11) 
2n \{l-x) + 2 J 

Forming forward matrix elements and setting the renormalization point /i = Q one gets 
from the RG equation fllCf ) the GLAP equation. Take into account Eq. ([TT1) one finds for 
x — > 1 

r \ q NS (y,Q 2 )-q NS (x,Q 2 " 



+2 r dy *~Kvw-<t"M)\ (12) 

Jx y — x I 



This equation is to be implemented by a initial condition that has, corresponding to Eq. 
(j|), the form q NS (x, Qq) oc (1 — x) v , where the value of v is approximately determined by 
the counting rule result. Using the ansatz q NS (x,Q 2 ) = A(Q,Q )(1 — x) u+ ^ Q ' Qo \ where 
£(Qo>Qo) — 0, it is straight forward to derive the solution WB 



q NS (x, g 2 ) oc - -ff^ n e M-**W°\l - xr^\ (13) 

where 

«nn\ or f Q2 dta s (t) AC F ( lnQ 2 /A 2 \ 

t{Q,Qo)=2C F J^ T — = — ^{^qTJ^), A, = 11 - (2/3W, (14) 

is a function increasing with Q. Thus, with rising Q the (1 — x) power behaviour increases 
and the prefactor decreases. 



4 Factorization of the distribution function 

In this section we present a more rigorous treatment that allows the perturbative computation 
of the large xsj behaviour of the parton distribution functions and also provides a smooth 
inclusive exclusive connection. As we saw in Section 2, at XBj — » 1, the parton model 
suggests that the momentum of the struck quark is far off shell. Thus, using the light-cone 
quantization and the representation of the wave function by a superposition of fock states, 
a factorization can be performed based on the iteration of the equation of motion || |] 

q(x, Q 2 ) ~ f 1 du I dv0*(u, \ 2 )T H (x, u, v, A 2 , Q 2 )<p(v, A 2 ), (15) 
Jo Jo 

3 A similar treatment is used for the factorization of exclusive large momentum transfer processes. 
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where the factorization scale is given by A 2 ~ m 2 /(l — x). Here the hard amplitude 
T H (x, u, v, A 2 , Q 2 ) includes the LC operator O a (x,fi = Q) and provides, in principal (see 
below), the power behaviour in (1 — x). 

For simplicity, in the following we restrict ourselves to the parton distribution functions 
of the meson. Since the large momentum arises from the kinematics of the partons inside 
the meson, the leading order approximation of Th(x,u,v, A 2 , Q 2 ) is determined by the two 
gluon exchanges (see Fig. 2). Note that Th is not dependent on the reference momentum 




to P V-U to 




P1 JU^l P2 




p to 



P-Pi P-k P-P2 



Fig. 1: This graph yields to the leading (1 — Fig. 2: Parton distribution function of the meson 
x) 3 power behaviour of the nucleon structure at large x. The crossed circle symbolizes the LC 
function for x — ► 1. operator insertion. 

Qo and that the Q-dependence will be induced by the renormalization of the LC operator. 
Furthermore, the meson distribution amplitude 0(w, A 2 ) satisfies an evolution equation which 
controls the logarithmical corrections in A 2 [HI: 



(u, A 2 ) =f>*(l- u)uCT{2u - 1) (|^|^ 2 



2^ \ t2/A> 

(16) 



k=0 



where 7° = Cp (3 + ( - fc+1 2 ^ fc+2 ^ ) — 4 YmZ\ 7) > C^ 2 are Gegenbauer polynomials, are nonper- 
turbative coefficients, and A 2 , = m 2 /(l — Xq) is an appropriate reference momentum square. 

Now we study the large x behaviour of the parton distribution functions for the scalar 
0L-theory. First we compute the hard amplitude T H by including the tree approximation of 
the LC operator vertex, which is given by 8{x — k + /P + ). Using the far off shell behaviour of 
the struck partons [see Eq. (f|)] one finds corresponding to Fig. 2 (substitute the gluon lines 
by scalar lines) the following Feynman integral representation, 

rr I \ f 1 A f 1 '" 1 A f 1 -* 1 -* 2 , l-X X -X 2 -X Zxt 

T H (x,u,v) ~ dx\ dx 2 dx 3 — _ 2 d(x - x\u - x 2 v - x 3 ), 

.xx^— - + x 2 ^—^ ~ M h j ~ [1 - (xi + x 2 )x 3 }m 2 , (17) 

where the mass of the spectator partons was set equal to the free particle mass m. The 
restriction of the integration region by the 5-function comes from the operator vertex and 
provides three contributions, 



Th(x, u,v) = 9 (u, v < x) T# + {# (u < x < v) + {u «-> u}| + 9 (x < u, v) T x 

T^(x,u,v) = (1 -x) (j— ~) (jZ ~) h( x > u > v )> 

7g(x,u,v) = (l-x)(- — -f 2 (x,u,v) + - — -f 2 (x,u,vj) , 

VI — u v — u J 

Th(x,u,v) = (1 - x)f 3 (x,u,v), 



3 

Hi 
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where fi, i = 1,2,3, and /2 are smooth functions, which can be more or less replaced by 
constants. 

The convolution (|1^) provides the x — > 1 behaviour; but because of the complicated 
structure of Th(x,u,v) in Eq. ([18]) the result depends on the end point behaviour of the 
distribution function <p(u). In the physical case, where the distribution amplitude cj>{u) 
vanishes at the end points, the leading term comes from the integration region u, v < x ~ 1, 
which provides^ 



q(x) = du dv(f)*(u, \ 2 )Tjj{x, w, v )4>{v , A 2 ) ~ (1 — x) c 



u 



for x -> 1. (19) 



It is also interesting to investigate the unphysical case, in which the distribution amplitude 
is constant or divergent at the end points, i.e., 4>{u) ~ (1 — u)~ e for u — > 1, e > 0. Then one 
gets for 1/2 > e R the behaviour 



q(x) 



:i-x) 



[x r ~ 



(l-a;)^. 



(20) 



Obviously, the obtained power behaviour is different from that given in Eq. ([19]). 

To take into account the evolution in Q we determine 0(x, fi 2 /k 2 , a s (/i)) for x — » 1 
from the RGE and include it in the computation of T#. For x — > 1 the integral kernel 
is completely determined by the anomalous dimension j^(a s ) of the field 0, P(x; a s ) = 
(a s / (27r))7?5(l — x) + 0(const., a 2 s ). Thus, the RGE is simplified to 



d 



0(x,/i /k } a s (/j,)) 



a s (jj? 



-^0(x^ 2 /k 2 ,a s (^)). 



d f i 2 v_ "" '"• 7 "^ r '" 2tt 
The solution of this first order (partial) differential equation can be found in a text book: 

0(x,fi 2 /k 2 ,a s (fi)) = 0(x,l,a s (k 2 )) 



(21) 



\\n{k 2 /A 2 ) J 



Six — 



k+ 



0{a s {k 2 )) 



/ ln(/x 2 /A 
\hx{k 2 /h? 



2\\ 27S//30 



(22) 



Here 0(x, 1, a s (k 2 )) represents the initial condition, which, assuming a s (k 2 ) is small enough, 
can be expanded perturbatively. 

Including the solution in the computation of Tjj and setting the renormalization 
point \i = Q provides for a physical distribution amplitude: 



g(x,Q s 



ln(g 2 /A 2 



27S//30 



x-t-1 



ln(m 2 /((l - x)A 2 



1 0(m, m 2 / (1 — x)) 



l-u 



(23) 



The comparison with the solution of the GLAP equation shows that the reference momentum 
square Ql = m 2 / (1 — x). 

4 Because of the space-time dimension six the obtained power behaviour is 3 and it agrees with the counting 
rule result. 

5 This restriction is necessary to avoid divergencies caused by the convolution. 
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5 Summary 



Using both the definition of the parton distribution functions as forward matrix elements 
of light-cone operators and the parton model describing the hadronic bound state we have 
factorized (for Q 2 > m 2 /(l — XBj) and large m 2 /(l — Xsj)) this distribution function as a 
convolution of a hard amplitude and a distribution amplitude. Here the factorization scale 
A ~ m 2 / (1 — XBj) originates from the kinematics of the parton inside the hadron as well as 
the assumption that the spectator mass does not vanish. Furthermore, we have computed 
the hard amplitude in the scalar toy theory for the distribution function of the "meson" 
and found that the resulting (1 — xsj) behaviour is governed by the end point behaviour of 
the distribution amplitude. For a physical distribution amplitude that vanishes at the end 
points, the obtained (1 — xbj) power behaviour agrees with the counting rule result. In the 
other case one finds a disagreement. 

Moreover, to obtain the evolution in Q we have included the renormalized LC operator 
in the computation of the hard amplitude. Since the initial condition of the renormalization 
group equation can be computed perturbatively (this is not possible for the initial condition 
of the GLAP equation) we have shown that the reference momentum square Qq is given by 
the xsj dependent factorization scale, i.e., by m 2 /(l — XBj)- 

Note that the use of this treatment for QCD is straight forward. It will be interesting 
to investigate the question of the (1 — xsj) power behaviour on the dependence of the 
distribution amplitude and the effects of higher twist terms. It is to be expected that the 
founded value Ql ~ m 2 /(l — Xbj) will be the same in QCD, so that the power behaviour of 
the valence quark distribution functions will be changed [see Eq. flT3|)]. 
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